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We derive a partial-wave (Mie) expansion of the axial 
force exerted on a transparent sphere by a laser beam focused 
through a high numerical aperture objective. The results hold 
throughout the range of interest for practical applications. 
The ray optics limit is shown to follow from the Mie expan- 
sion by size averaging. Numerical plots show large deviations 
from ray optics near the focal region and oscillatory behavior 
(explained in terms of a simple interferometer picture) of the 
force as a function of the size parameter. Available experi- 
mental data favor the present model over previous ones. 
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Optical tweezers are single-beam laser traps for neu- 
tral particles that have a wide range of applications 
in physics and biology Q. Dielectric microspheres are 
trapped and employed as handles in most of the quan- 
titative applications. The gradient trapping force is ap- 
plied by bringing the laser beam to a diffraction limited 
focal spot through a large numerical aperture microscope 
objective. 

Typical size parameters [3 = ka (a = microsphere ra- 
dius, k = laser wavenumber) range in order of magnitude 
from values < 1 to a few times 10 1 . A theory of the trap- 
ping force based on geometrical optics (GO) M should 
not work in this range. Other proposals (cf. pF), based 
on Mie theory, have employed unrealistic near-paraxial 
models for the transverse laser beam structure near the 
focus, incompatible with its large angular aperture. 

We take for the incident beam before the objective, 
propagating along the positive z axis, the usual Gaussian 
(TEM)oo transverse laser mode profile, with beam waist 
wq at the input aperture, where kwo 3> 1. We employ the 
Richards and Wolf j| representation for the correspond- 
ing strongly focused beam beyond the objective, with a 
large opening angle 9q (no paraxial assumption), taking 
due account of the Abbe sine condition. This should be 
a more realistic representation. 

The microsphere, with real refractive index (we ne- 
glect absorption), is immersed in a homogeneous medium 
with refractive index n±. We consider here the simplest 
situation, in which the sphere center is aligned with the 
laser beam axis, so that we evaluate the axial trapping 
force. With origin at the sphere center, we denote by 
r = — qz the focal point position. The fraction A of total 
beam power that enters the lens aperture is 



A = 1 - exp(-27 2 sin 2 



(1) 



where 7 is the ratio of the objective focal length to the 
beam waist wq. 

By axial symmetry, the trapping force in this situa- 
tion is independent of input beam polarization: we take 
circular polarization. The electric field of the strongly 
focused beam (we omit the time factor exp(— iuit)) has 
the Debye-type 0] integral representation 



E (r) =E Q 



2 7T 



d0 sin 9\J cos 9 exp (— 7 2 sin 2 1 



x exp [ik • (r + qz)] e(0, cj>), 



(2) 



where k — |k(0, (f>)\ — n.iw/c, e{6,<j>) = x! + iy', and 
the unit vectors x' and y' are obtained from x and y, 
respectively, by rotation with Euler angles a — <j), [3 = 
9, 7 = —(f). The factor V cos 9 arises from the Abbe sine 
condition. 

For each plane wave exp(ik • r) in the superposi- 
tion (^), the corresponding scattered field is given by 
the well-known Mie partial- wave series in terms of 
the Mie coefficients ai,bi, that are functions of the size 
parameter (3 and the relative refractive index n = nijn\. 
By substitution into (||), we obtain the total scattered 
field E s (r). 

The trapping force is found by replacing the total field 
E = Eo + E s (likewise for B) into the Maxwell stress 
tensor and integrating over the surface of the sphere. The 
resulting axial force F is proportional to the focused laser 
beam power P, 



F=( ni /c)PQ, 



(3) 



where Q is the (dimensionless) axial trapping effi- 
ciency 0. 

We denote by Q e the contribution from terms in EqE s 
and BqB s (that also give rise to the extinction efficiency) 
and by Q s the remaining terms, so that Q — Q e + Q s . 
We find 
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where G\ and Gi are multipole coefficients for the fo- 
cused beam, 



Gi = d6 sin 9V cos 9 exp(— 7 2 sin 2 9) exp(iS cos ( 
Jo 



x d[J6), 



Gt 



-idGi/dS, 



(6) 
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with (5 = fcg. In (||), are the matrix elements of 

finite rotations 0, that can be expressed as 



d[ 1 (6) = [pi(cos9)+U(cos9)}/(2l + l) 



(8) 



in terms of the Mie angular functions || pi and ti. The 
results (f|) and (^), apart from converging beam effects, 
have the same structure as the radiation pressure effi- 
ciency ||, with which they are closely related, as will be 
seen below. 

In the Rayleigh limit, (3 <C 1, Q is dominated by the 
electric dipole Mie term a\, and the trapping force (||) 
becomes F = (a/2)VE 2 , where a is the static polariz- 
ability of the sphere Q|. In the opposite limit /3 ^> 1, 
the connection with geometrical optics is established by 
applying to (^) and (||) the following steps |§ @- (i) In 
(6), substitute pi and ti by their (non-uniform) asymp- 
totic expansions for large I, and approximate G; and G l 
by the method of stationary phase |8) . (ii) Compute the 
average < Q > over a size parameter range associated 
with a quasiperiod of the Mie coefficients. The result is 



< Q >go-- 
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In (g), 9 1 = arcsin(gsin#/a), 9 2 = arcsin(sin 9\ /n) are 
the angles of incidence and refraction (defined so as to 
be negative if q < 0) at the sphere surface associated with 
a component in the direction 9 of the focused beam (j^) . 
The corresponding Fresnel reflectivity for polarization j 
(||, _L) is rj. Eq. (^|) may also be derived in the frame- 
work of GO. Thus, the expression within curly brackets 
agrees with the GO result for the force exerted by each 
component ray as first obtained in H. The remaining 



pre-factors in (^) , not accounted for previously, represent 
the intensity distribution of the focused beam as implied 
by the sine condition and the transverse profile of the 
laser beam at the input aperture of the objective. 

In Fig. 1, Q is plotted as a function of q/a, the center 
offset from the focus in units of the sphere radius p0[ . 
The numerical values chosen correspond to the experi- 



ment of Ref. (Tl| 
78°, which by~([[ 
represents the GO result 
sent the exact Mie results 



ni = 1.33, n 2 = 1.57, A = 0.85, O = 
yield 7 2 = 0.99. The dotted curve 
[). The other curves repre- 
1) and (|5t) for two different 



/3 values, corresponding to microsphere radii employed 
in [jllj: 1.42/im (dashed) and 2.16/mi (solid), respectively 
(3 = 18.8 and (3 = 28.4. The qualitative behavior of the 
GO curve has been explained in terms of competition 
between radiation pressure (scattering force) and gradi- 
ent force. However, the GO result for the maximum back- 
ward trapping efficiency Q m is smaller (by a factor of the 
order of 2) |l^] than the values obtained in Ref. This 
is in line with the discrepancy between experimental and 
theoretical values noted in Ref. |l3|| . 

The Mie theory provides values for Q m (0.088 for (3 = 
18.8 and 0.086 for [3 = 28.4 Q) below the GO result 
Q m — 0.095. The position at which the backward force is 
maximum lies beyond the corresponding GO value q/a = 
1.01. The stiffness decreases as this point is approached 
from the focus, contradicting the GO prediction and in 
agreement with an experiment reported in Ref. Jl3[ . 

GO is also a poor approximation near the geometrical 
focus, as expected. In fact, Fig. 1 shows that the exact 
values deviate substantially from GO near q = 0. The 
stable equilibrium position shows large positive as well 
as negative offsets from GO (further discussed below), 
and the linear Hooke's law range around the equilibrium 
position is narrower than predicted by GO. Because of 
the axial focusing effect || , the (nonuniform) asymptotic 
approximations to the Mie angular functions employed in 
the derivation of (|]) break down at q = 0, although (|]) 
is continuous at this point, yielding, with 9\ = 9 2 = 0, 



< Q >go (q = 0) 
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< cos 9 > , 
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where r is the Fresnel reflectivity for normal incidence 
and < cos 9 > denotes an average over the intensity dis- 
tribution of the focused beam (^). Since the incident 
rays are either backscattered or undeviated in this ap- 
proximation, ( [To| ) represents pure radiation pressure in 
GO. 

The region around q = deserves special treatment, 
in view of its relevance to the evaluation of trap axial 
stiffness. For ^> 1, the above discussion and the lo- 
calization principle imply that the main contributions to 
([|) and (^) should arise from partial waves with I <§; /3, 
so that we apply Hankel's asymptotic expansions to the 
spherical Bessel functions in the Mie coefficients. The 
results are independent of I, and the summations over 



2 



multipole coefficients can then be carried out, resulting 
in 



Q(q = 0) 



8rsin 2 A/2 
1 + r 2 — 2r cos A 



< cosv >, 



(11) 



where A = ^niioaj c. This expression corresponds to the 
radiation pressure efficiency (twice the reflectivity) of an 
infinite set of parallel-plate interferometers (width 2a, re- 
fractive index ri2, so that A is the round-trip phase), each 
one oriented at an angle 9 with respect to the axis, tra- 
versed at normal incidence by the corresponding beam 
angular component. The GO result ( jl(i| ) follows from jll| ) 
by taking an incoherent average. Since n — 1 is small, 
we have r <C 1, so that the interferometer reflectivity is 
nearly sinusoidal. 

In Fig. 2, for the same parameters as in Fig. 1, we 
plot Q at q = as a function of j3. The Mie curve (full 
line) displays the expected near-sinusoidal oscillation as 
(3 increases, approaching the interferometer behavior JTl| ) 
(shown in dotted line). The GO value (|l0|) (dashed line) 
is approached in the average sense. The two points corre- 
sponding to the P values employed in Fig. 1 are shown by 
circles. Since the radiation pressure at f3 = 28.4 is above 
the GO value, the Mie value for the equilibrium position 
q cq is larger than the GO result, in agreement with Fig. 1 
(the opposite applies at (3 = 18.8). The values for <7 eq are 
found by numerically solving the equation Q(q cq ) = 0. 
In the limit /? 3> 1, they are vanishingly small at [3 val- 
ues that are minima of Q(q = 0). For (3 ~ 5, q eq /a as a 
function of j3 oscillates in phase with the oscillations of 
Q(q = 0), around the GO value (g eq /a)GO = 0.217, and 
with amplitude of the order of 0.17. 

The trap axial stiffness is given by 



n x P ( dQ 



(12) 



Within GO, k decreases as 1/(3. This follows from scaling: 
Qgo depends on q only through q/a. Hence, dQco/dq = 
Q'Go(<l/ a )/ a > yielding 



niP 



(13) 



Again for the parameters of Ref. [[llj (power P — 3mW), 
we plot in Fig. 3 the Mie values of n (solid line) [jl5[ , 
the GO result kqo = (18//3)(pN//xm) calculated from 
Eq. (|ll]) (dotted line) and the experimental data points 
from Ref. pl| , with the respective error bars. We also 
show (dot-dashed line) the values predicted by the elec- 
trostatic model recently suggested by Tlusty et al. p6| . 
As could be expected, their approach may be applied 
only in the low- frequency (Rayleigh) limit, where it may 
be replaced by the simpler electric dipole approximation 
(neglecting the variation of the field over the sphere vol- 
ume) already discussed above in connection with (0) . In 



order to test the sensitivity of the results to the focused 
beam parameters, we also plot the Mie values for n cor- 
responding to a larger waist: 7 2 = 0.3 (dashed line). 

For (3 ~ 10, the Mie values for k oscillate around 
the GO curve with period A/3 = 7r/(2rt), like the force 
Q(q = 0) [cf. (JTT|)] and the equilibrium position. This cor- 
responds to a frequency interval Av = c/(4ri2a), which 
is in the THz range for spheres with radii of a few mi- 
crons. As shown in the inset in Fig. 3, where we plot 
(< k > —kgo)/kgo a s function of /3, the average of the 
Mie values, < n >, stands above the GO curve, but the 
relative difference decreases to zero as (3 increases beyond 
(3 « 55. For large sphere diameters, k may become very 
small over short (3 intervals. This may be of interest for 
applications to scanning force microscopy O]. 

In conclusion, by deriving an analytic Mie expansion 
for the axial trapping efficiency, based on a more real- 
istic model of optical tweezers, we are able to cover the 
range of interest for most applications. Furthermore, the 
connection with the correct GO limit (taking into ac- 
count the sine condition) has been derived from the Mie 
expansion by size averaging. The behavior near the fo- 
cus has been obtained and interpreted in terms of an in- 
terferometer model, which also accounts for equilibrium 
position and trap axial stiffness oscillations. These oscil- 
lations should be accessible to experiments by scanning 
the laser beam frequency. Most experimental data points 
lie above the GO values, closer to the wave optics predic- 
tions computed from the Mie expansion employing only 
the experimentally given parameters. 
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Fig.1 : Maia-Neto & N 

FIG. 1. Normalized axial force versus position (in units of 
the sphere radius) . The dotted line is computed from ray op- 
tics theory, whereas the solid and dashed lines are calculated 
from the wave-optics theory with size parameters (3 = 28.4 
and j3 = 18.8, respectively. The vertical dashed lines mark 
the microsphere boundaries. 




fig. 2: Maia Neto & Nussenzveig 



FIG. 2. Normalized force at the geometrical focal point ver- 
sus size parameter /3 : exact (full line); interferometer model 
(dotted line) and geometrical optics (horizontal dashed line). 
The black circles indicate the two values of /3 used in Fig. 1. 
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FIG. 3. Axial stiffness k of the optical tweezer as a func- 
tion of f3. Solid, dotted and dot-dashed lines correspond to 
the (exact) wave-optics theory, geometrical optics, and elec- 
trostatic theory [16], respectively, and for a focal length to 
waist squared ratio r y 2 — 0.99. Also shown are the experimen- 
tal data points of Ref. [11], with corresponding error bars, and 
the exact values for 'y 2 = 0.3 (dashed line). In the inset, we 
plot the relative discrepancy between the average of the exact 
values and the geometrical optics results (for -y 2 = 0.99). 
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